Nonequilibrium perturbation theory of the spinless Falicov-Kimball model 
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We perform a perturbative analysis for the nonequilibrium Green functions of the spinless Falicov- 
Kimball model in the presence of an arbitrary external time-dependent but spatially uniform electric 
field. The conduction electron self-energy is found from a strictly truncated second-order pertur- 
bative expansion in the local electron-electron repulsion U. We examine the current at half-filling, 
and compare to both the semiclassical Boltzmann equation and exact numerical solutions for the 
contour-ordered Green functions from a transient-response formalism (in infinite dimensions) on the 
Kadanoff-Baym-Keldysh contour. We find a strictly truncated perturbation theory in the two-time 
formalism cannot reach the long-time limit of the steady state; instead it illustrates pathological 
behavior for times larger than approximately 2/U. 
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I. INTRODUCTION 

The theoretical description of the physical proper- 
ties of strongly correlated electrons in an external time- 
dependent electric field is an important problem in con- 
densed matter physics. One reason is the potential for 
using strongly correlated materials in electronics, due to 
their tunability. These materials have unusual and in- 
teresting properties, which can be modified by slightly 
varying physical parameters, like temperature, pressure, 
or carrier concentration. For example, the conductance 
of nanoelectronic devices can be controlled by tuning 
the carrier concentration and applying an electric field 1 . 
Moreover, because of the small size of nanoelectronic de- 
vices, an electric potential of the order of 1 V produces 
electric fields on the order of 10 5 — 10 6 V/cm or higher. 
Therefore, it is important to understand the response of 
such systems to a strong electric field. Recent progress 
in nanotechnology has resulted in many experimental re- 
sults on strongly correlated electron nanostructures in 
time-dependent electric fields which must also be un- 
derstood. One such system is the quantum dot, which 
can be described by a Kondo impurity attached to two 
leads, through two tunnel junctions^. Other examples 
of important experiments on strongly correlated systems 
in time-dependent electric fields include an unusually 
strong change in the optical transmission in the transition 
metal oxide^ S^CuOa. and the dielectric breakdown of a 
Mott insulator that occurs in the quasi-one-dimensional 
cuprates S^CuOa, and SrCu02, when a strong electric 
field is applied 5 ". Generally speaking, since most electron 
devices have nonlinear current-voltage characteristics, it 
is important to understand how electron correlations will 
modify this effect. 

Due to the complexity of these problems, there is a 
dearth of theoretical results on different models avail- 
able. Similar to equilibrium, the simplest cases to an- 
alyze are the one-dimensional case and the approxima- 
tion of infinite dimensions, where dynamical mean-field 
theory (DMFT) can be applied. The zero dimensional 
problem (quantum dot) in the case of a Kondo impu- 



rity coupled to two leads, or sources of electrons, has 
also been studied by different groups. In particular, the 
case of a s ingl e-impurity Anderson model was examined 
in Refs. Ia l7l8l and a Bethe ansatz approach with open 
boundary conditions has now been developed for similar 
problems^. An attempt to describe optical transmission 
experiments 4 in SrCuOa was performed in the framework 
of the one-dimensional two-band Hubbard model, where 
the problem on a twelve site periodic ring was solved ex- 
actly. In Ref. 0, the spectral properties of the Hubbard 
model coupled to a periodically time varying chemical 
potential were studied to second order in the Coulomb 
repulsion with iterated perturbation theory. The possi- 
bility of the breakdown of a Mott insulator in the ID 
Hubbard model was discussed in Ref. [lfj by numeri- 
cally solving the time-dependent Schroedinger equation 
for the many-body wave function. The electrical and 
the spectral properties of the infinite-dimensional spin- 
less Falicov-Kimball model in an external homogeneous 
electric field have also been examinedii*i2*H. 

The spinless Falicov-Kimball model is one of the sim- 
plest models for strongly correlated electrons that ac- 
quires the main features of a correlated electron sys- 
tem. It consists of two kinds of electrons: conducting 
c-electrons and localized /-electrons. These two differ- 
ent electrons interact through an on-site Coulomb repul- 
sion. The model can be regarded as a simplified ver- 
sion of the Hubbard model, where the spin-up (c) elec- 
trons move in the background of the frozen spin-down 
(/) electrons. The Falicov-Kimball model was originally 
introduced as a model to describe valence-change and 
metal-insulator transitions 1 ^, and later was interpreted 
as a model for crystallization or binary alloy formation 15 . 
Much progress has been made on solving this model with 
DMFT in equilibrium, where essentially all properties of 
the conduction electrons in equilibrium are known (for a 
review, see Ref. GJ). 

We have already studied the nonequilibrium response 
of the conduction electrons in the Falicov-Kimball model 
to a uniform electric field turned on at a particu- 
lar moment of time with an exact numerical DMFT 
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formalismii*i2ii^. We solved a system of generalized 
nonequilibrium DMFT equations for the electron Green 
function and self-energy and found interesting effects 
caused by the electric field. In particular, we saw how 
the conduction-electron density of states (DOS) evolved 
from a Wannier-Stark ladder of delta functions at the 
Bloch frequencies to broadened peaks with maxima lo- 
cated away from the Bloch frequencies. Another interest- 
ing phenomenon we found was that the Bloch oscillations 
of the electric current can survive out to long times in the 
interacting case, when the amplitude of the electric field 
is large. These results were determined by numerically 
solving a system of the equations for the Green func- 
tion and self-energy defined on a complex time contour 
(see Fig. QJ, according to the Kadanoff-Baym-Keldysh 
nonequilibrium Green function formalism 1 -'" 1 -. The nu- 
merical results were obtained for a finite time interval 
because computer resources restrict the contour to be a 
finite length. Therefore, the steady state could not be 
rigorously reached. 



Our original motivation for this work, was to find a 
simple perturbation theory that could reach the steady 
state, at least for weak electron-electron interactions. 
Surprisingly, we found that a transient-based perturba- 
tion theory breaks down once the time becomes larger 
than 2/U, as described in detail below. Nevertheless, 
even though the perturbation theory cannot give a defi- 
nite answer about the steady state of the system, it does 
allow us to study the physical properties of the system 
in the transient regime, and for small U, we can extend 
the work farther than can be done with the exact numer- 
ical techniques. Therefore, it might help us understand 
the evolution of the system toward the steady state in 
the case when the electron correlations are weak. In this 
contribution, wc calculate the response of the Falicov- 
Kimball model to an external electric field by calculating 
the strictly truncated electron self-energy to second order 
in U. We focus on the time-dependence of the electric 
current and the density of states in the case of a constant 
electric field turned on at time t = 0. 



The paper is organized as follows: We define the prob- 
lem in Section II and introduce the nonequilibrium dy- 
namical mean-field theory formalism to solve it in Section 
III. In Section IV, the expressions for the nonequilibrium 
Green functions are derived. These functions are used 
to study the time-dependence of the electric current in 
Section V and the electron density of states in Section 
VI. Conclusions and discussions are presented in Section 
VII. 



II. FALICOV-KIMBALL MODEL 

The Hamiltonian for the spinless Falicov-Kimball 
model has the following form: 
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(i) 

The Falicov Kimball model has two kinds of spinless elec- 
trons: conduction c-electrons and localized /-electrons. 
The nearest-neighbor hopping matrix element ty = 
t*/2Vd is written in a form appropriate for the infinite- 
dimensional limit where d is the spatial dimension of the 
hypercubic lattice and i* is a rescaled hopping matrix 
element; note that our formal results hold in any di- 
mension, but all numerical calculations are performed in 
the infinite-dimensional limit. The summation is over 
nearest neighbors i and j. The on-site Coulomb repul- 
sion between the two kinds of the electrons is equal to 
U. The symbols [L and [if denote the chemical poten- 
tials of the conduction and the localized electrons, re- 
spectively. For our numerical work, we will consider the 
case of half-filling, where the particle densities of the con- 
duction and localized electrons are each equal to 0.5 and 
H = Hf = U/2, 

We study the response of the system to an external 
electromagnetic field E(r, t). In general, an external elec- 
tromagnetic field is expressed by a scalar potential tp(r, t) 
and by a vector potential A(r,<) in the following way: 



E(r,t) = -Vv»(r,t) 



1 dA(r, t) 
c di ' 



(2) 



For simplicity, we assume that the electric field is spa- 
tially uniform. In this case, it is convenient to choose 
the temporal or Hamiltonian gauge for the electric field: 
<p(r, t) = 0. Then, the electric field is described by a spa- 
tially homogeneous vector potential A(t), and the Hamil- 
tonian maintains its translational invariance, so it has a 
simple form in the momentum representation [see Eq. lfl| 
below]. This choice of the gauge allows one to avoid ad- 
ditional complications in the calculations caused by the 
inhomogeneous scalar potential. The assumption of a 
spatially homogeneous electric field is equivalent to ne- 
glecting all magnetic field effects produced by the vector 
potential [since B(r,i) = V x A(r,t)]. This approxima- 
tion is valid when the vector potential is smooth enough, 
that the magnetic field produced by A(r, t) can be ne- 
glected. When the electric field is described by a vector 
potential A(r, t) only, the electric field is coupled to the 
conduction electrons by means of the Peierls substitution 
for the hopping matrix 20 : 



Uj exp 



Uj exp 



He 



A(r,t)dr 



R, 



-A(t) • (R, - Rj). 



(3) 
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where the second equality follows for a spatially uniform 
electric field. 

The Peierls substituted Hamiltonian in Eq. (JJJ assumes 
a simple form in momentum space: 



h(a) = y: 



eA(t) 

he 



P,k,q 



(4) 



in Eq. the creation and annihilation operators are 
expressed in a momentum-space basis. The free electron 
energy spectra in Eq. (@J is 



k- 



eA(t) 

he 



(=1 



eA'(t) 

he 



(5) 



We shall study the case when the electric field (and 
the vector potential) lies along the elementary cell 
diagonal^: 



A(t) = 4(t)(l, !,...,!) 



(6) 



This form is convenient for calculations, since in this case 
the spectrum of the non-interacting electrons coupled to 
the electric field becomes 



e k 



eA(t) 
he 



en 



sm 



he J 

which depends on only two energy functions 
e(k) = -2i^cos(afc z ) 



s(k), 



(7) 



(8) 



III. NONEQUILIBRIUM FORMALISM 

The nonequilibrium properties of the system in Eq. Q 
can be studied by calculating the contour-ordered Green 
function: 

GHhM) = -*<T c c k tf(*i)4ff(*2)) (11) 
- l Tr{e-^- t — )f c exp[-i ^^/(t)]^^)^^)} 

= Xre -/3ff(-t max ) ' 

where the time integration is performed along "the inter- 
action" time contour depicted in Fig. ^ (see, for exam- 
ple, Ref. l23l) and the time ordering T c is with respect 
to times along the contour. The angle brackets indi- 
cate that one is taking the trace over states weighted 
by the equilibrium density matrix exp[— j3H (— t max )]/Z 
with Z = Trexp[— (3H(— t max )]; this is done because the 
system is initially prepared in equilibrium prior to the 
field being turned on. 



max 







max 



and 



e(k) = -2t^sin(afc'). 



(9) 



The former is the noninteracting electron bandstructure 
(in zero electric field), and the latter is the sum of the 
components of the noninteracting electron velocity mul- 
tiplied by —1. 

The diagonal electric field is especially convenient in 
the limit of infinite dimensions d — > 00, when the 
electron self-energies are local in space, or momentum- 
independent. Then the joint density of states for the two 
energy functions in Eqs. ijSJ and © can be directly found. 
It has a double Gaussian form for the infinite-dimensional 
hypercubic lattice^: 



P2(e,e) 



1 



• exp 



nt* 2 a d 

Below we use the units, where a 



e 



e 

^2 



(10) 



t* = 1. 



FIG. 1: Kadanoff-Baym-Keldysh time contour for the two- 
time contour-ordered Green function in nonequilibrium. The 
red line corresponds to the region of time where the electric 
field is nonzero (t > 0). 

The indices H and / on the right hand side of Eq. 1|11|) 
stand for the Heisenberg and Interaction representations 
for the electron operators. H{— t max ) corresponds to the 
Hamiltonian in the case of zero electric field. The expres- 
sion in Eq. (|llf> is a formal generalization of the equilib- 
rium formulas to the nonequilibrium case. Therefore, all 
the machinery of equilibrium quantum many-body the- 
ory, including Wick's theorem and the different relations 
between the Green functions, the self-energies and the 
vertex functions can be directly applied to the nonequi- 
librium case as long as we replace time-ordered objects by 
the contour-ordered objects along the contour of Fig. 
(see, for example, Ref. 1231 for an appropriate discussion). 

In this work, we use the two-branch Keldysh contour, 
which is appropriate when we directly solve the prob- 
lem on the lattice (rather than using DMFT techniques, 
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where the contour is necessarily more complicated) . This 
contour consists of two horizontal branches: one runs to 
the right and the other to the left. Such a time contour 
is needed to be able to calculate both time-ordered and 
so-called lesser or greater Green functions, because it al- 
lows us to use time ordering along the whole contour to 
determine those objects. Because of the time ordering, 
we can use Wick's theorem to evaluate different correla- 
tion functions for the perturbative expansion. In general, 
there are four different Green functions one defines by 
taking the time coordinates along the upper + or lower 
— branch. These are the time-ordered, lesser, greater, 
and anti-time-ordered Green functions. They satisfy 



= i(c 



^k(*l )*2 ) 

-i{Tc kH (ti)ct H (t 2 )}, (12) 

^k(*l )*2 ) 

i H (t2)c kH (h)), (13) 

^k(*l ! h ) 

-i<Ckjr(*i)4ff(*a)), (14) 
^k(*l > h ) 
= -i<Te kH (*i)4 ff (* 2 )>, (15) 

where the + and — signs indicate whether the real times 
t\ and t 2 he on the upper (positive time direction) or 
lower (negative time direction) branch of the time con- 
tour, the symbol T is ordinary time ordering and the 
symbol T is anti time ordering. The more familiar re- 
tarded and advanced Green functions satisfy 



Gl{h,h) 
G<{h,t 2 ) 

G k (h,t 2 ) 



(h,t 2 ) 



(ti,t 2 ) 



G^(ti,t 2 ) ~ G^(ti,t 2 ) 
-Gl(h,t 2 ) 
Gk(ti,t 2 ) - 
-Gl(h,t 2 ) + G<{t u t 2 ) 



+ G>(h,t 2 ), 
GKUM) 



and the so-called Keldysh Green function is 

G«(h,t 2 ) = G<(i 1 ,t 2 ) + G>(t 1 ,i 2 ) 
= Gl{t x ,t 2 ) + Gl{t u t 2 )- 



(16) 



(17) 



(18) 



It is sometimes convenient to introduce a two-component 
Fermi-field operator 



Pk(t) 



_ ( <*(*+) 
Ck(t-) 



4w 



(<i(t + ),ct(t-)) (19) 



with the components also defined on the upper and lower 
time branches of the time contour. In this case, the fol- 
lowing matrix Green function can be introduced: 



G k (ti,t 2 ) 



G k (t+,t+) G k (t+,t-) 
Gk(ti,t 2 ) G k (t 1 , t 2 ). 



(20) 



In order to transform the matrix in Eq. I|2U|) to 
a form more convenient for computation and of- 
ten used to study nonequilibrium problems, we make 



the combined Keldysh and Larkin-Ovchinnikov (LO) 
transformation 2325 : 



Gl(h,t 2 ) ^ Lf z Gl(t u t 2 )L\ 



(21) 



where L — (l/V2)(fo —vr y ), and fi,i — 0,x,y,z, are the 
Pauli matrices. Under such a transformation, the Green 
function in Eq. I|20[) acquires the following form: 



Gk(ti,t 2 



( G*(t u t 2 ) G«(h,t 2 
{ 



Gl{h,t 2 ) 1 (22) 



where the retarded, advanced and the Keldysh Green 
functions are also expressed as 



Gi 



l (h,t 2 ) = -i8(h - t 2 ) {{c kH (h),ci H (t 2 )}) , (23) 



Gt 



(h,t 2 ) = iOfa-h) ({c kH (ti),(i H (h)}) , (24) 



Gk(*l.*3) = -i([ckff(tl),4ff(*2) 



(25) 



The square (curly) brackets correspond to the commuta- 
tor (anti-commutator). In this representation, the time 
arguments of Green functions (|23|I - H25(I are real times. 

In the nonequilibrium case, all the information about 
the properties of a system is contained in two basic elec- 
tron Green functions. The most often used ones are the 
retarded function in Eq. H23J) and the lesser function in 
Eq. fT3|) . 

The retarded Green function contains information 
about the spectra of the system, and the lesser Green 
function describes the occupation of these states. These 
two functions form an independent Green function basis, 
and any other Green function can be expressed by means 
of these two functions; in the equilibrium case, only one 
Green function is independent (since the states are dis- 
tributed according to the Fermi-Dirac distribution). In 
particular, it is easy to find a useful relation, which fol- 
lows from the definitions in Eqs. f2*3^l . $2ty . f2*5|) and i fHfy : 



G£(h,t 2 ) = \ (G«(t u t 2 ) - G*(t u t 2 ) + G£(t u t 2 )) . 

The self-energy representation, which corresponds to 
Eq. J23 has the same form: 



p c( M _/E£(t 1( t a ) E£(ti,f 2 )\ 

^^-y ^(h,t 2 )J- 



(27) 



In this representation one can also apply the standard 
Wick rules to evaluate the different operator averages of 
the perturbative expansion. In fact, it is possible to show 
that the nonequilibrium and equilibrium field theories are 
formally equivalent if one introduces the time ordering 
along the contour instead of the usual "equilibrium" time 
ordering 2 *. 
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The Dyson equation which connects the contour- 
ordered Green function in Eq. 12-?l) and self-energy, re- 
mains valid: 



G' Oc v^c s~1C 
k L k b k 



(h,t 2 ), (28) 



where Gy(tx, t 2 ) is the nonequilibrium Green function in 
the case of no interactions, and the product of the three 
operators in the last term is a shorthand notation for an 
implicit matrix multiplication of the continuous matrix 
operators over their respective internal time coordinates 
(evaluated along the contour). Its components can be 
found analytically. For the conduction electrons, we have 



x exp 



-i I dte (k - eA(t)) 

Jto 



x exp 



dte (k-eA(t)) 



(29) 



(30) 



G™{hM) = i{2/[e(k)-/x(°)]-}lexpM )(t 1 -i 2 )] 



x exp 



dte (k-eA(t)) 



(31) 



where /[e(k) - ^°)] = 1/{1 + exp[/3{e(k) - /j( )}]} is 
the Fermi-Dirac distribution function for free electrons. 
The symbol //°' is the noninteracting chemical potential. 
While the /-electron Green functions are: 



Ff°(ti,t 2 ) = -te(t 1 -t 2 y^ }it ^ t2 \ 



i J k°(*i,*2) = *0(t 2 -ii)e i "/ ) <**-'») 



F^°(t 1) t 2 )=i(2n f -l)e 1 ^ 



(32) 



(33) 



(34) 



Here, up is the noninteracting chemical potential for the 
/-electrons and n/ is the concentration of /-electrons. 
As expected, the free /-electron Green functions are 
momentum-independent. Note that at half-filling, we 

have /i^ ) = nP — and rif = 1/2, so that the non- 
interacting Green functions simplify, and in particular 
F* = 0. 

The Dyson equation in Eq. (|28|) is equivalent to the 
following system of three equations: 

G£(ii,i 2 ) = G£°(*i,*a) + [G^G^(h,t 2 ), (35) 



G£(h,t 2 ) = G£\t u t 2 ) + [G^G£](h,t 2 ), (36) 

G«( tl ,t 2 ) = [l + G^G^il + ^G^h) 

+ [G^Gt]{hM)- (37) 



This is derived by simply writing down the Dyson equa- 
tion for every matrix component. The 21 component is 
trivial, since every 21 matrix component is equal to zero. 
We omit the internal time integrals (over the real-time 
axis) in Eqs. (|35[) - l|37l) for sake of brevity. 

In general, it is difficult to solve this system of equa- 
tions, especially in the nonequilibrium case when elec- 
tron interactions are present. However, for the Falicov- 
Kimball model the electron self-energy is momentum- 
independent through second order in the interaction U, 
which simplifies the analysis. In the next Section, wc 
shall present the perturbative solution of Eq. (|28[1 . or 
equivalently Eqs. I|35l) - (|37[1 . for the Green function. 



IV. PERTURBATION THEORY 

We will investigate a non-self-consistent perturbative 
expansion that is strictly truncated to second order in U. 
We perform the expansion directly for the lattice Hamil- 
tonian, which is worked out (for the equilibrium case) in 
the Appendix. As mentioned above, the first and second- 
order self-energies are local because the /-electron Green 
function is local in all dimensions. The expression for the 
nonequilibrium self-energy can be obtained from the cor- 
responding expression for the equilibrium time-ordered 
self-energy by making the Langreth substitution^ for 
the corresponding Green functions, which tells us to re- 
place all integrals over the real time axis by integrals over 
the contour, and to replace all time-ordered objects by 
contour-ordered objects. The equilibrium time-ordered 
self-energy has the following expression when strictly 
truncated to second order in U (see Appendix A) 

Eim(*i,ta) = S lm [S(t! - t 2 )(Un f - M + // 0) ) 

+ U 2 G^(t 1 ,t 2 )FH(t 2 ,t 1 )F[l(t u t 2 )] (38) 

where I and m denote lattice sites and we have suppressed 
the time-ordered superscript on all quantities. Using the 
fact that the product of the /-electron Green functions is 
nf(l—nf), yields the final expression for the second-order 
equilibrium self-energy 



KiihM) = S lm U 2 n f (l - n^GHhM). 



(39) 



Note that this expression cannot immediately be used 
to find the equilibrium retarded self-energy. In order to 
find that, we can determine the time-ordered self-energy 
on the imaginary time axis, Fourier transform to Mat- 
subara frequencies and then perform an analytic contin- 
uation to the real axis to find the retarded self-energy. 
An alternate method, which we will adopt here, is to 
find the equilibrium lesser self-energy from the nonequi- 
librium formalism (evaluated in the equilibrium limit). 
Then the equilibrium retarded self-energy follows by tak- 
ing the difference of the equilibrium time-ordered and 
lesser self-energies. 

In the nonequilibrium case, we use the Langreth rule^ 
to replace the integrals over real time by integrals over 
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the contour. Then by examining time arguments on each 
of the two branches, and following the same strategy as 
in the appendix, we find 

£ T(2) (ii,i 2 ) = U 2 G T0 (t% , t2)F T0 (fa , t-2 ) F T0 (t 2 , * i()40) 

^ >{2) {fa,t 2 ) = U 2 G >0 (fa,t 2 )F> (fa,t 2 )F< (t 2 ,t&2) 
X f M( tl ,h) = U 2 G* (fa,t 2 )F™(fa,t 2 )F™(t 2 ,fa(}l3) 

Each pair of products of /-electron Green functions is 
equal to n/(l — rif), so after performing the Keldysh and 
the LO transformations, we get 



/ £-R(2) ^K(2) \ 

( o E A( 2 ) ) (*ii*2) = Un f {l-n f ) 



(44) 



/ G R0 G K0 \ 

X V G A0 j ^■ 1 ' 2 ^ 



where all the Green functions are local. For instance, 

G fl0 (t 1 ,t 2 ) = l^G«°(i 1 ,t 2 ). (45) 



This approach also allows us to determine the equilibrium 
retarded self-energy by simply evaluating the results in 
the absence of an electric field. This is an example of 
how one uses the nonequilibrium technique to complete 
an equilibrium analytic continuation. 

Now we will focus on the half-filled case [where rif(l — 
rif) = 1/4, /i^ ) = /i^ ' = 0, and the first-order self- 
energy vanishes]. The formal expressions for the Green 
functions can be found from the appropriate Dyson equa- 
tions [Eqs. lj35|) - (j!T7|l ]. Here, we write down the results 
only for the half-filled case: 

G^ 2 \t u t 2 ) = G™( tl ,t 2 ) + [GfE^Gf }(fa,t 2 ), 



(46) 



Gt (2 \h,t 2 ) = Gf(fa,t 2 ) + [Gf E^Gf ](i 1; i 2 ), 

(47) 

G?< 2 >(tx,t a ) = Gr{tuh) + [G^ 2) G™ (48) 
+ G^ A ^G A0 + G£° Z K WG£°](t u t 2 ) 7 

where the second-order self-energies X^ 2 ) and the free 
Green functions G° are given by the half-filled case in 
Eqs. (gU) and (O-JSIJ, respectively. (In the general 
case, we need to also include the first-order self-energy 
and its iterated second-order contribution.) 

The free c-electron local Green functions (in the 
presence of a spatially uniform electric field), which 
are necessary to calculate the self-energy (|44|l . can be 
found by summing Eqs. H29fl - I|31[) over momentum [as in 



Eq. g5j] 21 : 

G m (t u t 2 ) = -i9(t 1 -t 2 y^ 0)iti - t2) 



x exp 



dt cos (eA(t)) 



1-2 



dt sin (e A (t)) 



12 



(49) 



G A0 (h,t 2 ) = ie{t 2 -h)e^ 0) ^-^ 
{ 1 



x exp 



ti \ 2 

cit cos (eA(t)) 



ii 

ti x 2' 

dt sin (eA(i)) 



12 



(50) 



G K0 (h,t 2 ) = ie^° } ^-^ 



x exp ■ 



dt sin (eA(t)) 



x / dep(e)[2/(e- M (°))-l] 



x exp \ — ie / dt cos (eA(t)) f , (51) 
Jt 2 



where 



P(e) = —j= e " 

\/lT 



(52) 



is the free electron density of states on the hypercubic 
lattice in infinite dimensions (the generalization to finite 
dimensions is obvious). At half-filling we have fi^ = 
in Eqs. 

Thus, we have obtained the formal expressions for 
the second-order nonequilibrium retarded, advanced and 
Keldysh Green functions l|35|) - (|37|l . with the free Green 
functions defined by Eqs. $Z%-^T$, gHJl-lEIl> and the 
self-energies in Eq. I|44|l . In the following Section, we 
shall use these solutions to study the time-dependence of 
the electric current in the case when an external time- 
dependent electric field [Eq. J()J)] is present. 

Before proceeding with the nonequilibrium solutions 
we would like to study some of the equilibrium properties 
of the retarded self-energy E. We check the perturbative 
result for the self-energy in Eq. (|44|l versus the exact 
DMFT results in equilibrium and estimate the values of 
U above which the second order equilibrium perturbation 
theory fails. It is natural to suppose that the second 
order nonequilibrium perturbation theory results are also 
inaccurate when U is that large. 

Using the results in Eq. (|44|) . we find the expression 
for the second order self-energy has the following form in 
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frequency space (at half- filling) : 

El 

4 

El 

4 



— iirp(uj) 



(53) 



where the symbol P denotes the principle value of the 
integral, and p(e) is defined in Eq. (|52l) . In order to 
estimate the values of U, for which the second order per- 
turbation theory gives good results in equilibrium, we 
compare the imaginary part of the perturbative retarded 
self-energy 



ImS R (w) 



(54) 



with the exact self-energy calculated by the DMFT ap- 
proach. The frequency dependence of — ImE fl for differ- 
ent values of U is presented in Fig. As follows from 
this figure, the second order perturbation theory gives 
reasonable results up to U ~ 0.5. Hence, we expect that 
the result in Eq. I|44|) will also be accurate when U < 0.5 
for the noncquilibrium case. 



0.5 - 



hh 0.25 



- U=0.25 
,' \ - 0.50 
0.75 




FIG. 2: Frequency dependence of the imaginary part of the 
second order retarded self-energy for different values of the 
Coulomb repulsion U (solid lines). The dashed lines are the 
corresponding exact dynamical mean-field theory results. 



V. THE ELECTRIC CURRENT 

A. Semiclassical Boltzmann equation 
approximation for the current 

Before analyzing the time dependence of the cur- 
rent in the quantum case, we present the approxima- 
tion for the current that follows from a semiclassical 
Boltzmann equation approach. The Boltzmann equation 
for the nonequilibrium quasiparticle distribution function 

/ non (M) 



-iG^(t,t) becomes 



<9/ non (k,<) 
Ft 



- e E(t).V k / non (k,t) 



1 



(M)-/(k)], 

(55) 



in the case of a spatially uniform time-dependent elec- 
tric field. Here, r is the scattering time, which typically 
is proportional to l/U 2 in the weak correlated limit of 
the Falicov-Kimball model (because it is inversely pro- 
portional to the imaginary part of the self-energy at zero 
frequency) and the r.h.s. of the equation is the collision 
term. Note that in the quantum case, the quantum Boltz- 
mann equation is more complicated because the quantum 
excitations depend on two independent times, not one. 
Hence, it isn't obvious that we should use the quantum 
relaxation time in the semiclassical Boltzmann equation; 
indeed, we often find that by fitting to a somewhat longer 
relaxation time, we can improve the fit to the Boltzmann 
equation rather dramatically, but we present the results 
using a fixed relaxation time here. Of course, our numeri- 
cal nonequilibrium results presented below are equivalent 
to exactly solving the full quantum Boltzmann equation, 
but our method of solution is different. 

The boundary condition for the semiclassical Boltz- 
mann equation is: 



r on (k,f = o) = /(k) 



i 



exp[/5(e(k)- /i)] + l' 



(56) 



We will be considering only the case of half-filling, where 
we can set the chemical potential to zero. The solution 
of Eqs. Ij55(l and (|56|) when the electric field is pointing 
in the diagonal direction E(t) = E(l, 1, l)0(t) has the 
following form: 

/ non (k, t) = ~ J die- (t -^ /T f (^ + e J t dt'-E(t'^J . 

(57) 

The total current can now be determined from the dis- 
tribution function 

j{t) = -~ X>(k) J* dte-^' T f (k + ej* dt'E(t')J 

(58) 

Shifting the momentum via k — > k — e J dt'E(t') and 
then integrating over the Brillouin zone gives the follow- 
ing analytical expression for the current: 



m 



eEr 



dep(e)ef(e) 



1 + e 2 £ 2 r 2 
1 - (cos(eM) - e£>rsin(e£^))e~ t/7 



(59) 



Analysis of this solution shows that the current is a 
strongly oscillating function of time for t -C r (when E 
is large) which then approaches the steady-state value 



■steady 



eEr 
l + e 2 £ 2 T 2j0 ' 



where 



3o 



dep(e)ef(e), 



(60) 



(61) 
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for t > r (see Fig. [3] for an example of the solutions in 
the infinite-dimensional limit). It is interesting that the 
amplitude of the steady-state current is proportional to 
E in the case of a weak field (the linear-response regime), 
and then becomes proportional to l/E when the field is 
strong (eE 1/t). In this nonlinear regime, the ampli- 
tude of the current goes to zero as the field increases. As 
it will be shown in the following Subsection, the exact 
quantum solution of the problem results in qualitatively 
different behavior from the Boltzmann case for the cur- 
rent as a function of time. 




the expressions for the retarded, advanced and Keldysh 
Green functions on the r.h.s. of Eq. i|26|) are given by 
Eqs. 146(1 - 1|48|1 . Since the electric field lies along the lat- 
tice diagonal, all its components are equal to each other, 
and the magnitude of the electric current becomes: 



(64) 



for any component a The momentum summation in 
Eq. Hti3|) can be performed by using the definitions of 
energy functions in Eqs. (JSJ) and @ and the expression 
for joint density of states in Eq. JTUJ: 



x [ecos (eA a (t)) — esin (eA a (t))] 

x [G* (t, t) - G* (t, t) + G*(t, t)] . (65) 



The sum of the last two terms in Eq. (|65J) gives zero, 
since 



FIG. 3: Boltzmann equation solution for the time-dependence 
of the electric current for different values of the electric field 
in the limit of infinite dimensions at r = 1, /3 = 10. Here and 
in Figs. we set the prefactor e/\fd in the results for the 
electric current [see Eqs. 15911 and 1681 1 to be equal to 1. 



G*SiM + GiSxM = i<{ck(ti),<£(*2)» 

= i5(h-t 2 ) (66) 



B. Perturbative expansion for the current in the 
quantum case 



is momentum-independent and the velocity is an odd 
function of momentum. 

Therefore, 



In the quantum case, the electrical current can be 
found by calculating the expectation value of the elec- 
trical current operator. This operator is formed from the 
sum over momentum of the product of the electric charge 
times the velocity vector at k times the number operator 
for electrons in state k. Hence, the time-dependence of 
the expectation value of the electric current components 



= e E -^rj ( k - eA W) <4(*M*)>> ( 62 ) 



which becomes 



(t) = -i-j= J2 sin ~ eA » (*)) G k (*. «). (63) 



m = 



de J dep2(e,e) [e cos (eA a (t)) 



2Vd. 

esin(eA a (t))} [G^(t,t) 

^ J — oo J — oo 

<R0r, , \r<R0/, , \r*K0i 



x {G^{tM)GZ{tiM)G^{t 2 ,t) 
+ G™{tM)Gt°{tx,t 2 )G^{t2,t) 
+ G^(t,t 1 )Gf o c (ti,t 2 )Gf i (h,t))], 



(67) 



after using the definition of the lesser Green function. 
This Green function can be found from Eq. (|26|l . where 



where we restricted to the half-filling case, used Eq. I|44|l. 
and neglected the term that is fourth order in U. Inte- 
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FIG. 4: The PT electric current as a function of time for 
E — 1, f3 = 10 and different values of U. Note the unphysical 
increase in the current at long times. 



gration over e in Eq. (|65|l yields 



Vd 
eU 2 



dep(e)ef(e) sin(eA(t)) 



4Vd J- 



dti 



esin(eA(i)) 



exp 



dt 2 — S(t 2 ,ii)cos(eA(t)) 



{C 2 (t 2 ,t 1 ) + 2S 2 (t 2 ,t 1 )} 



x J dep{e)[2f{e) - 1] exp[-ieC(t 2l h)] 

ieU 2 /"* /"* 

dti / <ft 2 [5(t 2 ,ti)cos(e J 4(t)) 



16-v/d 

- C(t 2 ,ti)sin(eA(t)) 
x exp -i{C* 2 (t 2 ,ii) + 25 2 (t 2 ,ii)} 
x / dep(e)[2/ (e) - 1] ex P [ieC(t 2 , t x )] 



(68) 



In order to simplify the expression, we have introduced 
the functions 

ta 



(69) 



C(i 2 ,*i) = / di cos(eA(t)) 
S{t 2 ,h) = / disin(eA(t)). 



In our numerical work, we consider the current in the case 
of a constant electric field E turned on at time t = 0, i.e. 
A(f) = —Et6{t). In this case, we find: 



C(t 2 ,h) = 



(-ti)0(-t 2 ) [ta - *i] 
sin(e-Ei 2 ) 



0(-ti)O(t 2 ) 

e(h)0(-t 2 ) 



-ti 



fa - 



sin(ei?ti) 



S(t2,*l 



0(h)6(t 2 )— [sm(eEt 2 ) - sm(eEh)} (70) 
-0(-ti)0(t 2 )-^[l-cos(e£7t 2 )] 

e(ti)5(*2)-= [oos(eMi) - cos(e£i 2 )](.71) 
eh 



Equations l|68|l . (|70|l and l|71|l determine the time de- 
pendence of the current when a constant electric field is 
turned on at time t = 0. It is difficult to find an analytic 
expression for j(t), except for the simplest case of U = 0. 
In this case: 



j(t) = j Q sin (eEt) 



(72) 



(the so-called Bloch oscillation 2 ^) , where the amplitude 
of the oscillations jo is given by Eq. (|61|l and the Bloch 
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FIG. 5: The PT electric current as a function of time for 
E — 0.5, P — 10 and different values of U (red lines). The 
exact DMFT calculation results and the Boltzmann equation 
solution are presented by black and green lines, respectively. 



FIG. 6: The PT electric current as a function of time for 
E — 1.0,/3 — 10 and different values of U (red lines). The 
exact DMFT results and the Boltzmann equation solution 
are presented by black and green lines, respectively. 



frequency is^l u>b = eE. In the limit of infinite dimen- 
sions, the current amplitude jo can also be written as: 



]a = \Tdl dep{€) 



de 



(73) 



Comparison of this result with the general expression 
for the current given in Eq. 1)68(1 yields the following for- 
mal expression for the current in a strictly truncated per- 
turbation theory expansion: 



j{t) = j sin (eEt) + U 2 j 2 (t). 



(74) 



after integrating by parts. 



Thus, the electric current is a superposition of an oscil- 
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lating part and some other part, whose time-dependence 
will be determined below. 
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FIG. 7: The PT electric current as a function of time for 
E = 2.0,(3 = 10 and different values of U (red lines). The 
exact DMFT calculation results and the Boltzmann equation 
solution are presented by black and green lines, respectively. 
Note that the perturbative values for the current become so 
large that they are not shown for larger times. 



Numerical results for the time-dependence of the elec- 
tric current calculated from Eq. I|fci8fl ar e presented in 
Figs. 0H3 As follows from these figures, the current os- 



cillates for all time. This is the main difference from 
the Boltzmann equation result (Fig. |3J) , where the cur- 
rent rapidly reaches a constant steady-state value. As 
it will be shown below, the second order perturbation 
theory cannot describe the steady-state regime correctly. 
However, the perturbation theory should give reasonable 
results for times smaller than ~ 2/J7, which we consider 
below. As depicted in Figure^ the current depends only 
weakly on U, except for the case U ~ 0.5. This is be- 
cause the current is dominated by the Bloch term j [see 
Eq. I)74p]. As time increases, the amplitude of the os- 
cillations decrease. This decrease is proportional to U 2 , 
in agreement with Eq. iJBSJ. As t increases further, the 
current amplitude goes through zero (at a time we call 
tpeit), and then starts to increase dramatically. This is 
an unphysical result, so we assume that the perturbation 
theory is not valid for times larger than t pcrt , i.e. the time 
when the perturbation theory breaks down. This time is 
always smaller than the time needed to see the steady 
state develop, so this perturbation theory is not capa- 
ble of producing the steady-state response. The period 
of the oscillations remains essentially equal to the Bloch 
oscillation period, but it is not well defined for damped 
oscillations. 

The validity of the second-order noncquilibrium per- 
turbation theory can also be checked by comparing the 
results for the current with exact numerical results calcu- 
lated by the nonequilibrium DMFT method 11 ! 12 ; 13 . We 
present the corresponding results for cases of different 
values of U and E in Figs. |SHZ| As shown in these Fig- 
ures, the perturbation theory results are close to the ex- 
act results only for times smaller than ~ 2/U. Note that 
although the perturbation theory does show an increase 
in the current for intermediate times, it does not prop- 
erly show the quantum beats in the current, which occur, 
with a beat period on the order of 1/U, even though the 
shape of the curve is qualitatively similar. These beats 
occur only for large fields (E > 2 here). 

In Figures |SHZI we also present the results for the 
current in the case of the Boltzmann equation solution 
Eq. I|59|) . where the scattering time is fixed at the quan- 
tum Boltzmann equation value r — /o(/x)/(47r|ImS(o; = 
0)|) (see, for example Ref. |3 this result is the trans- 
port relaxation time). Substitution of Eq. 1(31}) into this 
expression yields: 



= l/(w 2 U 2 ). 



(75) 



The Boltzmann equation approach shows a rapid ap- 
proach to the steady state, much faster than what is 
seen in the exact numerical results or in the perturbation 
theory. Although the Boltzmann equation is accurate at 
small times, it is clear that the quantum system has much 
richer behavior than what is predicted by the semiclas- 
sical approach, at least in cases where the electric field 
is large. Note that the functional form of the semiclassi- 
cal Boltzmann equation can fit the exact results for the 
current much better if we adjust the relaxation time r to 
yield the best fit of the data rather than use Eq. (|75|l , but 
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the relaxation time then becomes just a fitting parameter 
and is not derived from a microscopic model. 

To complete the analysis of the time-dependence of the 
electric current, we present an analytical expression for 
the electric current in the limit of strong electric fields 
E -> oo. Substitution of Eqs. (ZD) and (HI into Eq. itfiSl 
gives the following result for the electric current in the 
limit of large electric fields: 



m 



-^=J dep(e)ef(e) 
l-U 2 B([3)-^-t 2 



sm(eEt), 



(76) 



VI. LOCAL DENSITY OF STATES 

The time-dependent density of states can be calculated 
from the local retarded Green function in the following 
way: 



A(u,T) = IuiG r (uj,T), 

IT 



(78) 



where 



where G r (uj,T) is the Fourier transform of the retarded 
Green function G R {ti,t<z) with respect to the relative 
time coordinate r = t± — t2, and T is the average time 
coordinate T = (t\ + t2)/2. The perturbation theory 
retarded Green function is determined in Eq. (|46|l . In- 
tegration over the Brillouin zone and substitution into 
Eq. CBJ yields: 



B(/3) 



dep(e)ef{e) 



dx I dyye v 



x / dep(e)f(e)sm(2sy). (77) 



Numerical calculations show that this is a positive de- 
creasing function of temperature: 0.25 < B{(3) < 0.5. 
As follows from Eq. (|76() . the electron-electron corre- 
lations lead to a decrease of the current amplitude at 
short times, since the terms proportional to U 2 in the 
square brackets in Eq. H76fl have a minus sign in front 
of it. This correction is not large at short times, since 
we consider the case U <C 1. However, from the nu- 
merical calculations, we find that the term in Eq. I|76|) 
which is proportional to J7 2 i 2 /4, is dominant in the 
case of longer times. In fact, the current is an oscil- 
lating function of time. The oscillation amplitude, de- 
creases with time as 1 — U 2 B{(3) — U 2 t 2 /A. At the time 
t pert = {2/U)s/i - U 2 B{(3) ~ 2/U the sign in front of 
the amplitude changes, and the current oscillates out of 
phase with the U — case. It is important that the 
period of the oscillations is equal to the period of the 
Bloch oscillations in the noninteracting case. At longer 
times, the term proportional to t 2 is the most important 
and the time dependence is j(t) ~ t 2 s'm(eEt) at times 
longer than ~ 2/U. These results are in a qualitative 
agreement with the numerical results for the perturba- 
tion theory presented in Figs. 0][7| However, it is plain 
to see that the a growing amplitude for the current as 
time increases for t > 2/U is an unphysical result. 

Unfortunately, it is impossible to find an analytical 
expression for the current in the limits of intermediate 
and small fields. However, our numerical analysis shows 
that the current qualitatively has the dependence given in 
Eq. H76fl . In particular, the steady state is never reached 
before the perturbation theory develops pathological be- 
havior. 



A(u,T) 



x exp 



dr cos(wt) 



1 



__(C 2 (T + t/2,T-t/2) 



+ S 2 (T + t/,T-t/2)) 



jj2 pT+t/2 pt 3 

x <! 1 / dt 3 dt 4 



T-t/2 



T—t/2 



x exp 
x exp 



"2 (C 2 {UM) + S 2 {UM)) 



(C(T + T/2,T-r/2)C(t 3 ,t 4 ) 



+ S{T + t/2,T -T/2)S(hM))]}, (79) 



where C(t 3 ,t4) and S(tz,t±) are defined in Eqs. I|70|l and 

CD- 

Numerical calculations of the time dependence of the 
density of states show that the PT correction is small for 
times less than t pclt ~ 2/U. Similar to the caseSi U = 0, 
the density of states develops peaks at the Bloch frequen- 
cies. These peaks grow as time increases, essentially due 
to the zeroth-ordcr DOS. However, these peaks do not 
split into two peaks as time increases for times less than 
~ 2/U, except for the zero frequency peak. This result is 
different from the exact DMFT calculations 12 at U = 0.5. 
Generically, one expects the DOS to broaden as scatter- 
ing is added, and to split, by an energy of the order of U 
in the steady state due to the interactions; the numeri- 
cal calculations indicate that this is indeed what actually 
happens even though the truncated perturbation theory 
does not explicitly show it. 

It is interesting that one can also find analytical expres- 
sions for the zeroth and the first two spectral moments, 
which are time independent (the first and second-order 
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moments are written for the half-filling case): 

J duiA(u,t) = 1, 

J dujujA(uj,t) = 0, 

J dcJw 2 A(tJ,t) = ^ + ^. (80) 

These results coincide with the exact results obtained in 
Ref. E3 at half-filling. 

VII. CONCLUSIONS 

We have studied the response of the conduction elec- 
trons in the spinless Falicov-Kimball model to an exter- 
nal electric field by using a strictly truncated perturba- 
tive expansion to second order in U . We derived ex- 
plicit expressions for the retarded, advanced and Keldysh 
Green functions and used them to calculate the time- 
dependence of the electric current and the density of 
states. We examined the case when a constant electric 
field is turned on at a particular moment of time. We 
found that Bloch oscillations of the current are present 
up to times at least as long as t ~ 2/U, but with a de- 
creasing amplitude. In the perturbative approach, the 
oscillations are around a zero average current and the 
PT breaks down for longer times, so we cannot reach the 
steady state. This result is different from that of the 
semiclassical Boltzmann equation approach. There, the 
Bloch oscillations rapidly decay with time and the system 
approaches the steady state with j = const, as time goes 
to infinity. Increasing the Coulomb repulsion reduces the 
amplitude of the oscillations more rapidly. The pertur- 
bation theory and the Boltzmann equation results for the 
electric current are close for short times. 

The perturbation theory results for the Green func- 
tions can be used to test the accuracy of the numerical 
nonequilibrium DMFT calculationsii*i2ii^ in the case of 
weak correlations, when the system is not in the insulat- 
ing regime. We also found excellent agreement for short 
times, but the PT becomes ill behaved at long times. 

Our results can be generalized to more complicated 
cases. In particular, cases with different time-dependence 
for the field (pulses, periodic fields etc.) can be consid- 
ered. Also higher order terms can be taken into account 
in the PT, or one can make the PT self-consistent. Fi- 
nally, we can examine other strongly correlated models 
as well, where we expect similar behavior. 

The most important result of this work is that it shows 
how difficult it is to accurately determine the steady- 
state response in a perturbative approach. At the very 
least, one needs to perform a self-consistent perturbation 
theory to be able to properly determine the long-time 
behavior, but, since we know that equilibrium PT is most 
inaccurate at low frequencies, it is natural to conclude 
that nonequilibrium PT is most inaccurate at long times. 



Hence, our work shows explicitly how challenging it is to 
try to determine the steady states of quantum systems 
unless one can formulate a steady-state theory from the 
start. Work along those lines is in progress for the exact 
solution via DMFT. 



Acknowledgments 

V. T. would like to thank Natan Andrei for a valu- 
able discussion. We would like to acknowledge support 
by the National Science Foundation under grant num- 
ber DMR-0210717 and by the Office of Naval Research 
under grant number N00014-05-1-0078. Supercomputer 
time was provided by the DOD HPCMO at the ASC and 
ERDC centers. We would like to thank the Referee of our 
paper—, which appeared in the proceedings of the Work- 
shop on Progress in Nonequilibrium Physics III (Kiel, 
Germany, August, 2005), for suggesting that we exam- 
ine the Boltzmann equation solution for the current and 
compare it to solutions found with other methods. 



APPENDIX A: EXPANSION FOR THE 
TIME-ORDERED EQUILIBRIUM SELF-ENERGY 
TO SECOND ORDER IN U 

The time-ordered self-energy E/ m (ii,i2) can be ex- 
panded by examining an order-by-order solution of the 
Dyson equation (we suppress the T superscript in these 
formulas) 

G^t') = G%{t,t') 

+ Y.J dtl f dt 2 G il (t,t 1 )^i m (t 1 ,t 2 )G mj (t2,t'lAl) 

lm 

= G%{t,t') 

+ £ /#i / dt 2 G^(t,t 1 )El2(tut 2 )G° mj (t 2 ,t') 

lm 

+ E/*i / *2GS(*.*i)S I ^(ti,ta)GS» i (<2,* / ) 

lm 

+ J2 J dtl J dt2 J dt3 J dt ±Gti(t,ti) 

Imnr 

x sW(t 1 ,t 2 )G^„fe,i 3 )Sl 1 J , ) (t3,t4)G° J (t4,t'), (A2) 

where the second equation explicitly shows the system- 
atic expansion to order U 2 in terms of the first-order E^ 1 ) 
and second-order E( 2 ) perturbative self-energies. The 
self-energies are found by a straightforward perturbative 
expansion of the evolution operator, when expressed in 
the interaction picture: 

iG l3 (t,t') = E^-J- / dtl ~ f dt » 
i/=0 v - J J 

x (f[H u (t 1 )...H u (t„)c t (t)c](t')}) con (M) 
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where the interacting part of the Hamiltonian is 

i i 

+ Uj^flM* (A4) 

i 

and the statistical average is taken with respect to the 
free Hamiltonian H = -J2uj) <w4 c j ~ A* (0) J2i c i c i ~ 

[if J2i fl fi- The subscript "con" indicates that we take 
only the connected contractions generated by Wick's the- 
orem, which must connect the conduction electron oper- 
ators for all time values in the matrix elements; the con- 
nected diagrams are included because the disconnected 
diagrams cancel from a similar expansion of the parti- 
tion function, which appears in the denominator of the 
matrix-element average. We have included the terms 
with the shift of the chemical potentials into Eq. HA4[) . 



since (p, — fi^) ~ ([if — M/ ) ~ U. Expanding the elec- 
tron Green function in Eq. (|A3|) up to second order in 
the interaction and using the definition of the self-energy 
in Eq. (|A2|) . one finds 

£ffi(ti,t2) = {Un f - f i + [i^)S lm S(t 1 -t 2 ), (A5) 

and 

£;™(*i,*2) = U 2 G° m (h - t 2 )Fi m {ti - t 2 )F® nl (t 2 - h), 

(A6) 

which depends on the difference of the times because 
we are in equilibrium. Since the /-electron Green func- 
tion is local, the second-order self-energy must have 
/ = 777. Furthermore, since all of the Green functions 
are time-ordered Green functions, the product of the two 
/-electron Green functions is equal to 77/(1 — nf). 
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